Superoscillation in speckle patterns 
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Waves are superoscillatory where their local phase gradient exceeds the maximum wavenumber in their Fourier 
spectrum. We consider the superoscillatory area fraction of random optical speckle patterns. This follows 
from the joint probability density function of intensity and phase gradient for isotropic gaussian random wave 
superpositions. Strikingly, this fraction is 1/3 when all the waves in the two-dimensional superposition have 
the same wavenumber. The fraction is 1/5 for a disk spectrum. Although these superoscillations are weak 
compared with optical fields with designed superoscillations, they are more stable on paraxial propagation. © 
2008 Optical Society of America 

OCIS codes: 260.3160 (interference), 030.6140 (speckle), 050.4865 (optical vortices) 



The spatial rate of change of a plane wave is determined 
by its wavevector, which can be defined as the gradient of 
its phase. Its direction and magnitude (wavenumber) are 
unambiguous. The phase gradient of more complicated 
wave fields can be used as a definition of local wavevec- 
tor, which may be a complicated function of position. 
Such is the case in superpositions of plane waves. The 
modulus of phase gradient is sometimes smaller than the 
superposition's maximum wavenumber, and sometimes 
bigger. This latter case has drawn much attention: at 
such places, the phase changes more rapidly than the 
constituent plane waves, hence the term 'superoscilla- 
tion' [1]. Superoscillatory waves, which locally vary much 
faster than their fastest Fourier component , have surpris- 
ing and counterintuitive properties, and have recently 
been studied in a variety of systems, particularly signal 
processing, quantum mechanics and optics [2-5]. Quan- 
tum mechanically, they fit into the general notion of weak 
measurements [6], and applications in optical imaging 
science have been suggested [7]. 

Our purpose here is to study some simple superoscil- 
latory aspects of two-dimensional random waves, that 
is, superpositions of plane waves whose direction (in a 
plane) and phase are independent and uniformly distrib- 
uted random variables. Such superpositions are a well- 
established model for speckle patterns - scalar waves, 
either optical [8] or acoustic [9], reflected or refracted 
from random rough surfaces. In this case, the two dimen- 
sions we consider are those of the plane transverse to the 
overall propagation, and the waves are superoscillatory 
in the sense that the transverse phase gradient is larger 
than the maximum transverse wavenumber. These su- 
perpositions are also used as a model for quantum wave- 
functions in two-dimensional chaotic enclosures ('chaotic 
billiards') [10, 11]; in particular, when the system is not 
time-reversal symmetric (either due to absorption, open 
channels or a magnetic field), the wavefunction is com- 
plex. 



A significant fraction of the area of a typical speckle 
pattern, whose transverse wave spectrum is band- 
limited, is superoscillatory. Specifically, a wave tp = 
pexp(ix), dependent on planar position r — (x,y), with 
intensity I — p 2 and phase gradient Vx> is 

superoscillatory where |V%| 2 — fc^ ax > 0, (1) 

where fc max is the maximum wavenumber in the trans- 
verse superposition spectrum. This definition of super- 
oscillation originated in a recent study of the relationship 
between waves and rays in structured refractive materi- 
als [12]. The purpose of this paper is to examine the areas 
where Eq. (1) is satisfied in random optical waves. 

Optical vortices - the nodes of wave fields, where the 
phase is undefined [14] - may be thought of as extremes 
of superoscillation, since the phase gradient diverges as 
I —> 0. Clearly, the vortices lie in superoscillatory re- 
gions. The configuration of vortices in random waves and 
speckle patterns has been much studied [8,15-17], and 
the present work may be thought of as generalizing this 
to all parts of the wave where Eq. ([IJ is satisfied. 

The phase gradient V% is also related to the current 
density J = Im^*V^ = IV x- Therefore J and V% are 
parallel, but have different lengths; as will be demon- 
strated, superoscillation tends to occur where / is small, 
so the fluctuations in the |V%| distribution are greater 
than those for J = \J\. 

The calculation of the superoscillatory fraction of a 
random wave superposition will use gaussian statistics, 
as usual in the study of speckle patterns [18], which ap- 
plies in the limit of infinitely many independent ran- 
dom plane waves. For isotropic random two-dimensional 
waves, the complex fields ip, d x ?p and d y ip have indepen- 
dent gaussian probability density functions, with vari- 
ances (|^| 2 ) = Iq (mean intensity), {\djij)\ 2 ) = 2I k 2} 
where j = x,y and k 2 is the normalized second moment 
of the power spectrum (which is circularly symmetric in 
fc-space), with well-defined fc max - A natural measure of 
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the correlation length of the random field is k 2 . 

We begin with the calculation of the joint probability 
of intensity and current P(I, J), from its Fourier trans- 
form, 



P(I,J) = 



ds / dHexp(isI + it ■ J) 



(2) 



(2tt) 3 

x (exp(-is|^| 2 + \t ■ (tf>*Vip - ipVip*))) 



where (•) denotes the gaussian average. This average 
is straightforward to calculate: in the exponent, the 
quadratic forms depending on s and t may be added 
to those from the gaussian probability density, yielding 
a complex six-dimensional quadratic form matrix with 
determinant (2(1 + il s) + I$k 2 \t\ 2 ) 2 / I^fcf. The gaus- 
sian average is the reciprocal square root of this deter- 
minant, divided by the square root of the product of 
the variances, (I^k^/A) 1 / 2 . The integrals in s and t may 
then be found using straightforward complex integration 
techniques, and then the (equidistributed) direction of J 
may be integrated. The final result is 



P(I,J) 



J 



exp 



~(2I + J 2 /Ih) 
zi 



(3) 



Integrating over J gives the well-known distribution for 
intensity, P{I) = exp(—I/I )/I , and integrating over I 
gives the previously derived probability density for J in 
two dimensions, P(J) = 2JK (y / 2/k 2 J/I )/I$k 2 (cf. [9] 
Eq. (84), [11] Eq. (18)), where K is a modified Bessel 
function. 

Since J = J|Vx|, the joint probability density function 
for / and |Vx| c & n be found by dividing by the jacobian 
determinant \d(I, J)/d(I, | VxD! = I, 



P(I,\Vx\) = I -S^e X p 
ttk 2 



-^(l + |Vx| 2 /2fc 2 )|- (-11 



This probability distribution green is one of the main 
results of this Letter, and is plotted in Fig. [l] It shows 
quantitatively a clear correlation between low intensi- 
ties and high phase gradients (and vice versa) in ran- 
dom waves, in line with common wisdom on superoscilla- 
tion [12]. Integrating over I gives the probability density 
function for |Vx|: 



4fc 2 |V X | 



(2fc 2 + |V X | : 



(5) 



This simple expression was previously derived in [8] 
(Eq. 4-190; also see [13] for a more extended discussion). 
The probability distribution of the phase gradient is un- 
bounded, with diverging variance. Since, for gaussian 
random waves, area averages are equivalent to ensemble 
averages [18], the area fraction / of the speckle pattern 
which is superoscillatory is 



/ = 



d|V X |P(|V X |). 



(6) 
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Fig. 1. (Color online) Contour plot of the joint proba- 
bility density function P(I, | VxD °f Eq. Q. It is clearly 
unbounded both in / and |Vx|, although high values of 
these arc anticorrelated. 



The numerical value of the superoscillatory area frac- 
tion / depends on the transverse spectrum, specifically 
the relationship between fc max and k 2 . Mathematically, 
the most natural spectrum to choose is monochromatic 
waves in the plane. This corresponds to random wave 
solutions of the Helmholtz equation V 2 ^ + k 2 ip = (V 2 
is the two-dimensional laplacian operator). For paraxial 
beams, this corresponds to nondiffracting speckle pat- 
terns: the spectrum lies on a ring of radius k in trans- 
verse Fourier space, and in this case k 2 = k 2 = fc 2 ^. 
Putting this into Eq. ^ gives / ring = 1/3 : a third of a 
random monochromatic two-dimensional wave superpo- 
sition is superoscillatory. A numerical realization of such 
a random wave field is shown in Fig. [2] Clearly, all of the 
vortices are in the superoscillatory region. 

It is instructive to compare the superoscillatory ar- 
eas with the area of 1/3-lowest intensity, as in Fig.[2](b). 
Although the areas are similar, the contour |V%| = k en- 
closes a slightly different region from the intensity con- 
tour, and their topologies are not the same. Superoscil- 
latory regions are therefore subtly different from low- 
intensity regions. 

It is interesting to observe that, on dividing by exp(i^), 
the real part of the Helmholtz equation can be rewritten: 



|Vx| 2 -fc 2 = 



V 2 p 

p 



(7) 



Therefore, for general monochromatic waves, sub- 
and superoscillation is governed the laplacian of the 
wave's real amplitude. In particular, the superoscillatory 
boundary contours are given by the nodal lines of V 2 /?. 
It is interesting to observe that the right-hand side of 
Eq. Q is the negative of the 'quantum potential' in the 
hydrodynamic interpretation of quantum mechanics [19] . 
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(a) (b) 

Fig. 2. (Color online) A random superposition of 100 
two-dimensional plane waves with the same wavenum- 
ber k. (a) Phase pattern (hues); (b) intensity pattern 
(grayscale). The white contour denotes the line |Vx| = k. 
The suboscillatory region, occupying statistically 2/3 of 
the area, is shaded with a dark filter in (a) . Several phase 
singularities can be seen in the supcroscillatory region of 
(a). The intensity contour (dashed cyan) enclosing the 
lowest 1/3 of the intensity pattern is also shown in (b). 
This is close to the white superoscillation contour, but 
the two are clearly different. The area plotted is (Air/k) 2 . 

The superoscillatory fraction for other band-limited 
isotropic spectra is easy to calculate. For instance, for 
a top-hat spectrum (disk spectrum) with equal weight- 
ing for all waves with wave number k < fc max , it is 
easy to show that k% = k 2 lllx /2 [18]. The superoscil- 
latory fraction here is then /disk = 1/5. Interpolat- 
ing between the two is the case of an annular spec- 
trum, where fc max (l ~ 8) < k < k max , where S is a 
scaled thickness: 8 — 1 is the disk spectrum, and the 
limit <5 — ^ is the monochromatic ring spectrum. For 
the annular spectrum, k 2 = &^ ax (2 — 25 + S 2 )/2, and 
/annuius = 1 — 4/(6 — 26 + S 2 ) , which smoothly interpo- 
lates between the two limiting cases. 

Of course, several familiar speckle spectra, such as 
a gaussian distribution of fc, are not band-limited, and 
so, strictly speaking, cannot be superoscillatory. Math- 
ematically, in the paraxial regime - appropriate for 
speckle patterns - all transverse wavenumbers, even in- 
finitely large ones, are infinitesimal compared to the 
z-componcnt of the full wavevector. The full three- 
dimensional wavenumber is technically infinite, and 
there are no true paraxial superoscillations. However, the 
transverse wave, considered as a superposition of two- 
dimensional waves, is superoscillatory as we describe. 

It is, of course, possible to study superoscillation in 
volumes of three-dimensional wave fields, not subject 
to the drawback of paraxiality (such as the fields rele- 
vant to [12]). A natural choice would be superpositions 
of isotropically random monochromatic waves in three 
dimensions, modelling, for example, (scalar) field modes 
of chaotic cavities [10,17]. These calculations, general- 
ized to D dimensions, appear in a follow-up paper [20]. 

The stability of transverse, highly superoscillatory 



fields on paraxial propagation has been previously con- 
sidered [3]: the superoscillations were found to propa- 
gate towards the far field, but ultimately were suppressed 
by suboscillatory regions. The naturally-occurring super- 
oscillations considered here have far smaller phase gra- 
dients than those in [3] , but their area fraction is statis- 
tically constant on propagation. This observation origi- 
nates from the fact that, within the paraxial approxima- 
tion, the power spectrum does not change on propaga- 
tion. In the extreme case of transversely monochromatic 
speckle patterns, the Fourier spectrum itself is invariant 
on propagation - the beam is diffraction-free - and so 
the superoscillatory regions are themselves invariant on 
propagation. 

We are grateful to Michael Berry, John Hannay, Kevin 
O'Holleran and Miles Padgett for discussions. MRD and 
JC are supported by the Royal Society of London. 
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